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ABSTRACT 

It has been suggested that hypernova remnants, with a substantial amount of energy in semi- 
relativistic ejecta, can accelerate intermediate mass or heavy nuclei to ultra-high energies and provide 
sufficient amount of energy in cosmic rays to account for the observed flux. We here calculate the 
expected energy spectrum and chemical composition of ultra-high energy cosmic rays from such semi- 
relativistic hypernovae. With a chemical composition equal to that of the hypernova ejecta and a flat 
or hard spectrum for cosmic rays at the sources, the spectrum and composition of the propagated 
cosmic rays observed at the Earth can be compatible with the measurements by the Pierre Auger 
Observatory. 

Subject headings: cosmic rays — gamma rays: bursts 



1. INTRODUCTION 

Ultra high energy cosmic rays (UHECRs) are the most 
energetic particles ever detected in the universe. There 
is a consensus that UHECRs with energies larger than 
10 19 eV are of extra-galactic origin. However, the com- 
position of ultra high energy cosmic rays at these en- 
ergies remains disputed. Recent observations of their 
maximum air shower elongations < X max > and fluctua- 
tions RMS(< X max >) by the Pierre Auger Observatory 
(PAO) suggest that UHECRs are progressively domi- 
nated by heavier nuclei at energies above a few EeV 
(Abraham et al. 2010), while the results from High 
Resolution Fly's Eye Experiment (HiRes) and Telescope 
Array (TA) are consistent with a pure proton compo- 
sition (Abbasi et al. 2004; Tsunesada et al. 2011). 
These claims, however, depend on the poorly-understood 
hadronic interaction models at such high energies. Nev- 
ertheless, the advantage of heavy composition for UHE- 
CRs is that heavy or intermediate-mass particles are 
much easier to be accelerated to UHE energies than pro- 
tons as the maximum particle energy is rigidity depen- 
dent (i.e. proportional to the nuclear charge Z of the 
particles). Heavy nuclei with energy above 100 EeV 
will suffer from strong photo-disintegration attenuation 
due to interactions with cosmic microwave background 
(CMB) and cosmic infrared background (CIB) photons, 
with an attenuation length shorter than ^100 Mpc (see 
e.g. Puget et al. 1976; Allard et al. 2006; Hooper et 
al. 2007; Allard 2009; Aloisio 2011; Taylor et al. 2011). 
Thus, sources producing UHECRs above 100 EeV must 
be within ~ 100 Mpc. Within this distance, there are 
few sources that are powerful enough to be able to ac- 
celerate particles to energies 10 20 eV. The candidates 
include local active galactic nuclei (AGNs) (e.g. Bier- 
mann & Strittmatter 1987; Berezinsky et al. 2006; Pe'er 
et al. 2009) and gamma-ray bursts (GRBs) (e.g. Wax- 
man 1995, 2004; Vietri 1995; Wick et al. 2004; Dermer 
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& Atoyan 2006; Murase & Nagataki 2006; Murase et al. 
2006). 

Extragalactic hypernova has been suggested to be a 
candidate source of UHECRs as well(Wang et al. 2007, 
2008; Liu et al. 2011). Hypernovae are a peculiar type of 
supernovae that are more bright and have larger explo- 
sion energies than typical ones (Paczyhski 1998; Soder- 
berg et al. 2006). Several hypernovae have been ob- 
served in the nearby universe to date, usually associated 
with sub-energetic GRBs. For example, SN 1998bw, as- 
sociated with GRB 980425, is the first-found peculiar 
supernova at a distance of 38 Mpc, with an isotropic- 
equivalent total kinetic energy of ~ 5 x 10 52 erg (Galama 
et al. 1998; Kulkarni et al. 1998; Iwamoto et al. 1998). 
Radio afterglow of this event showed that more than 
10 50 erg of kinetic energy was released in the form of 
a mildly relativistic ejecta. SN 2006aj, associated with 
GRB 060218, was detected at a distance of 140 Mpc 
(Campana et al. 2006; Mirabal et al. 2006; Soderberg 
et al. 2006; Mazzali et al. 2006). The strong thermal X- 
ray emission component in GRB 060218 is interpreted as 
arising from a semi-relativistic supernova shock break- 
out, in which a mildly relativistic ejecta with energy 
about 10 50 erg has been inferred (Campana et al. 2006; 
Soderberg et al. 2006). Radio observations of the recently 
discovered supernova SN2009bb, for which no associated 
GRB is found, also suggest that more than 10 50 erg is 
coupled into the mildly-relativistic ejecta(Soderberg et 
al. 2010). We will use the term semi-relativistic hy- 
pernovae to denote such supernovae harbouring mildly 
relativistic ejecta. Such semi-relativistic ejecta drives 
semi-relativistic shocks, which could accelerate parti- 
cles to ultra- high energies (Wang et al. 2007). The 
recently proposed engine-driven supernova as a source 
for UHECRs are identical to the hypernovae in nature 
(Ray & Chakraborti 2011), since both are powered by 
the central engine and both have much more energy 
in semi-relativistic ejecta. The estimated event rate of 
SN2009bb-like engine-driven supernovae is also compa- 
rable to the hypernova rate(Soderberg et al. 2010). The 
progenitors of hypernovae are believed to be Wolf-Rayet 
(WR) stars, which are stripped of their original H-rich 
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and sometimes even He-rich envelopes. The circum- 
stellar winds are therefore rich in intermediate mass ele- 
ments, such as C and O. Moreover, the hypernovae ejecta 
contain heavy elements such as Fe that are produced by 
nuclear fusion in the core of exploding WR stars, as well 
as elements heavier than Fe that are synthesized during 
the explosions. 

In this paper, we study the propagated energy spec- 
trum and chemical composition of UHECRs that origi- 
nate from hypernovae, and confront them with the mea- 
sured data by PAO. The rest of this paper is organized as 
follows. In §2, we first briefly review the hypernova sce- 
nario and describe the source composition of UHECRs 
produced by hypernovae. In §3, we calculate the energy 
spectrum and composition of UHECRs arriving at the 
Earth after propagation and then confront them with the 
Pierre Auger observations. Finally, we give our discus- 
sions and conclusions in §4. The effect of intergalactic 
magnetic fields is neglected in this paper. Throughout 
the paper, we use eV as the unit of particle energy and 
use c.g.s units for other quantities and denote by Q x the 
value of the quantity Q in units of 10 x . 

2. THE CHEMICAL COMPOSITION AND ENERGY 
SPECTRUM OF UHECRS AT THE SOURCE 

2.1. Hypernova model for UHECRs 

Wang et al. (2007, 2008) have suggested that hyper- 
nova remnants can accelerate heavy or intermediate mass 
nuclei to ultrahigh energies above 10 19 eV. Here for sim- 
plicity we only focus on the semi-relativistic part of the 
ejecta, since only this part is relevant to the acceleration 
of UHECRs. Particles are accelerated in the shock re- 
gion where the semi-relativistic ejecta is freely expanding 
before being decelerated by the swept-up circum-stcllar 
medium. The size of this free-expansion phase region for 
ejecta of a particular velocity /3 s hC and kinetic energy Ek 
is i?HN = E k v w /T 2 c 2 M ~2x 10 17 -E fei 5ir- 2 Mr5<V3cm, 
where T = (1 — (3 2 h )~ 1/>2 is the bulk Lorentz factor of 
the ejecta, M — 10 _5 M_5MQyr _1 is the wind mass 
loss rate (whose average value is 3 x 10 -5 Moyr -1 for 
WR stars), and v w = 10 3 v w 3kms _1 is the wind velocity 
(Willis 1991; Chevalier & Li 1999). We assume a frac- 
tion of eb of shock internal energy goes into the mag- 
netic field and hence the magnetic field energy density is 
B 2 /8ir w 2T 2 l3 2 e B p w c 2 , where p w = M/AirR^Vw is the 
mass density of the stellar wind at radius Run- Thus, by 
equating the adiabatic cooling time with the acceleration 
time, the maximum energy of accelerated particles is 
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where Z is the nuclear charge number of the particle and 
e is the charge of electrons. Note that, in the semi- 
relativistic hypernova scenario, the adiabatic cooling of 
UHE nuclei is much more efficient than other cooling 
processes, such as synchrotron cooling, photopion pro- 
duction and photo-disintegration processes (Wang et al. 
2008), so we do not consider them here. 

Hypernova can also provide sufficient energy in UHE- 
CRs(Wang et al. 2007; Liu et al. 2011). According to 



Katz et al. (2009), the energy production rate 3 in UHE- 
CRs per logarithmic energy interval inferred from the 
measured flux by PAO is ~ 10 43,5 ergMpc _3 yr _1 . Hy- 
pernova ejecta with T/3 ~ 0.5 can produce cosmic rays 
with energies ~ 10 20 eV and the kinetic energy in such 
velocity ejecta is typically ~ 10 51 erg (e.g. Wang et al. 
2007; Chakraborti et al. 2011). Assuming that half of 
the kinetic energy goes into cosmic rays and that 10% of 
this energy is distributed in every logarithmic energy in- 
terval, the required local event rate of hypernovae is then 
~ 600Gpc~ 3 yr _1 , which is consistent with the observed 
rate (e.g. Soderberg et al. 2006; Liang et al. 2007). 

2.2. Composition and spectrum at the sources 

It is usually assumed that shock acceleration leads to a 
power-law energy spectrum for particles. The differential 
number of particles accelerated by one hypernova can be 
described by 
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where the subscription i represents the species of the 
nuclei (e.g. z=H, He, and etc.), Ei, max and £', l min are re- 
spectively the maximum and minimum energy of cosmic 
rays injected into the intergalactic space by the hyper- 
nova, and Ui is the total energy of nuclei of species i 
released by one hypernova into the intergalactic space. 
According to Eq. 1, hypernova remnants can accelerate 
iron nuclei to energies about - 5 x 10 20 eV, so we fix 
£ maX; Fc = 10 20 5 (Z/26)eV in the following calculation. 
-Ei,min is the low energy cutoff in the spectrum, since 
nuclei with energy below Ei^ m { n « ZieBhost-Rhost can not 
escape from the confinement by the magnetic fields in the 
host galaxies, where -Bhost is the typical magnetic field 
of the host galaxy and JShost is the scale height of the 
host galaxy. Taking our own galaxy as an approxima- 
tion, we have E mi ^ Fe « 10 19 (Z/26)B_ 6 (^/0.3kpc)eV, 
where h z is the scale height of the Galaxy. For a sin- 
gle velocity ejecta , the power-law index is usually as- 
sumed to be p ~ 2 for both non-relativistic shocks and 
semi-relativistic shocks (Kirk & Schneider 1987). In the 
following calculation, we will assume the spectrum to be 
power-law with exponential cutoffs above -E max and be- 
low E min . 

Given that both Sj imax and Ei jla i n are rigidity depen- 
dent (i.e. £j imax(min ) oc Z^H.ma^min)), the total mass 
of nuclei of species i released into the intergalactic space 
by one hypernova is 



/■Si, max jju-. a . 

M,, CR = Ai / ^-dEi oc ^ 



(3) 



3 The value obtained by Katz et al. (2009) is based on pure- 
proton composition assumption. Since iron nuclei attenuation 
length around 10 20 eV is similar to that of proton, we adopt this 
value as the energy production rate in pure heavy nuclei scenario. 

4 It has been suggested that hypernova produces a kinetic energy 
distribution spreading over a range of velocity as E k oc (rf3)~ a 
(e.g. Soderberg et al. 2006). The final injected CR spectrum from 
the hypernova should be contributed by different velocity ejecta 
and the superimposed spectrum becomes softer(Wang et al. 2007). 
We here consider the highest energy cosmic rays, and for simplicity, 
just focus on the fastest (semi-relativistic) part of the ejecta and 
assume p ~ 2 in the following calculation. 
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Approximating Zi = Ai/2 (except for H), and assuming 
that the same proportions of particles of each species get 
accelerated (i.e. the value of M^cr/M, is independent 
of i), we have 

U H :U i :U J = 2M s :M i :M J . (i,j^H,i^j) (4) 

With the above equation, one can transform the abun- 
dances of nuclei of each species to the abundances of 
these nuclei at certain given energy, as described by 
Eq.(2). 

In the hypernova scenario, cosmic ray particles origi- 
nate from the circum-stellar wind material or hypernova 
ejecta material. The WR stellar wind contains abun- 
dant intermediate mass elements such as C and O, which 
dominate the chemical enrichment of interstellar medium 
(ISM) (e.g. Abbott 1982; Bieging 1990). The mass ra- 
tio of different elements for typical WR stellar wind is 
M Hc : M c : M : M x = 0.32 : 0.39 : 0.25 : 0.04 (Bieg- 
ing 1990), derived based on observation of emission line 
of these elements in WR stars and the stellar evolution 
model (see e.g. Willis 1982; Prantzos et al. 1986; van der 
Hucht et al. 1986, where X denotes the elements other 
than He, C and O. In the following calculation, we will 
use this composition as a representative case for the cos- 
mic rays from the stellar wind. 

On the other hand, the hypernova ejecta contains both 
intermediate mass elements and heavy elements. It con- 
sists of products of nuclear reaction in the interior of the 
WR star, such as C, O, Mg, Fe, as well as some heav- 
ier elements produced by the explosive nucleosynthesis. 
Chemical composition of the ejecta can be derived by 
modeling the spectra and light curves of hypernovae (e.g. 
Nakamura et al. 2001; Fryer et al. 2006) or through nu- 
merical simulation of SN explosion based on some spe- 
cific stellar model (e.g. Georgy et al. 2009). Best fit for 
the early spectra and light curve of SN 1998bw is ob- 
tained in the model named "CO138E50" by Nakamura 
et al. (2001), where the input explosion kinetic energy 
and ejecta mass are, respectively, Ek — 5 x 10 52 ergs and 
M e j = 10M Q . This model yields a chemical composition 5 
of M c : M : M Ne : M Mg : M Si : M s : M Ca : M Fo = 
0.006 : 0.71 : 0.037 : 0.034 : 0.083 : 0.041 : 0.007 : 0.09. 
In the following calculation, we will use this composition 
as a representative case for cosmic rays from the hyper- 
nova ejecta. In both scenarios, the amount of hydrogen 
is negligible, because the progenitor of Type Ib/c SN is 
stripped of H envelope, as indicated by lack of H line 
in the SN spectrum (Iwamoto et al. 1998; Nakamura et 
al. 2001). So during the explosion stage, there is little 
amount of H in both the wind and the SN ejecta. 

3. THE PROPAGATED SPECTRUM AND CHEMICAL 
COMPOSITION 

3.1. Propagation of UHE CR nuclei 

While UHECRs are propagating in the intergalactic 
space, there are three main attenuation processes due to 
interaction with CMB or CIB photons: the photodisintc- 
gration process, the Bethe-Heitler process and the pho- 
topion production process. The photodisintegration pro- 
cess causes nuclei to lose their nucleons and thus change 

5 We neglect the elements with mass less than 10 — 3 Mq and 
treat elements heavier than Fe (e.g. Ni, Co) as Fe because such 
elements would decay into Fe finally. 



their species, but does not reduce their Lorentz factor, 
while the latter two processes can reduce the Lorentz fac- 
tors of nuclei. Thus, we need to consider the evolution of 
the mass number A(t) and Lorentz factor jN(t) of a nu- 
cleus with time jointly. In the hypernovae scenario, given 
that E ma x/A < 10 19 eV, the corresponding maximum en- 
ergy of each nucleon is below the threshold energy for 
photopion production with CMB photons even at high 
redshifts. Although these nuclei can interact with more 
energetic CIB photons, the photopion energy loss rate 
with such photons is much lower than the energy loss 
via the Bethe-Heitler process by interacting with much 
denser CMB photons, so we neglect the photopion energy 
loss in the following calculation. In this work, we adopt 
the CIB model suggested by Razzaque et al. (2009) and 
Finke et al. (2010), and use the online data of the CIB 
intensity at different redshifts 6 . 

When an UHE nucleus is propagating in the intergalac- 
tic space, its Lorentz factor evolves with time as 



dt 



= "/N(t)H(z)+j N . BH (t) 



(5) 



where r fN(t)H(z) represents the adiabatic energy losses 
due to cosmological expansion and 7at,bh(£) is energy 
loss rate due to the Bethe-Heitler process. Here H (z) = 

-ffo\/^m(l + z) 3 + 17a is the Hubble constant at time 
t = t{z), where H Q = 71kms" 1 Mpc" 1 , fl m = 0.27 and 
£1a = 0.73. In an isotropic photon background, the 
Bethe-Heitler energy loss rate of a nucleus of the Lorentz 
factor ji is given by 
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where 7 Pl Bii(i) is the Bethe-Heitler energy loss rate for 
a proton of the same Lorentz factor 7jv, Z is the nuclear 
charge of the nucleus, ctbh and /bh are the cross section 
and the fraction of energy loss in one interaction for the 
Bcth-Heitler process (Chodorowski et al. 1992), and eth 
is the threshold energy for the interaction. Here e 7 is 
the photon energy in the rest frame of the nucleus while 



e is the photon energy in the lab frame. n 7 (e, 



the 



number density of background (CMB+CIB) photons of 
energy e at redshift z. 

The nucleus suffers from loss of nucleons by the photo- 
disintegration process, so the mass number evolves as 



dA „ , 
where Ra is the reaction rate, which is 



(7) 
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where o-^m.A is the total photo-disintegration cross sec- 
tion for a nucleus of mass number A. The photo- 
disintegration cross section is dominated by the giant 

6 http://www.phy.ohiou.edu/~finkc/EBL/indcx.html 
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dipole resonance (GDR) up to 30MeV, with the thresh- 
old energy between about lOMeV and 20MeV for all nu- 
clei (in the nucleus rest frame) . The cross section in this 
energy range can be modeled by a Gaussian form (e.g. 
Puget et al. 1976; Stecker & Salamon 1999; Hooper et al. 
2008) or a Lorentzian form (e.g. Khan et al. 2005; An- 
chordoqui et al. 2007; Hooper et al. 2007). From 30MeV 
to 150MeV, the quasi-deuteron process becomes dom- 
inant and the cross section can be approximated as a 
plateau (e.g. Puget et al. 1976; Stecker & Salamon 1999; 
Ahlcrs & Taylor 2010). In this work, we use the tabu- 
lated cross section data generated by the code TALYS 
7 (Goriely et al. 2008) which considered all the individual 
nucleon emission channels for nuclei with A > 5. For 
nuclei with A < 5, we adopt the Gaussian form near 
the threshold along with a plateau at higher energies to 
describe the total cross section. 

We use the Runge-Kutta method to solve Eq. 5 and 
Eq. 7 jointly so that we can trace the evolution history 
of the Lorentz factor jn(z] z s , E Sl A s ) and mass number 
A(z; z s , E s , A s ) of a specific nucleus injected from the 
source at redshift z — z s , where E s and A s are the initial 
energy and mass number of the nucleus. The evolution of 
the mass number is assumed to develop along the Puget- 
Stecker-Brcdckamp chain (see e.g. Puget et al. 1976; 
Stecker & Salamon 1999; Ahlers & Taylor 2010). We 
also record the redshifts, Lorentz factors and number of 
secondary protons 8 that are disintegrated from the par- 
ent nuclei so that the contribution of secondary protons 
to cosmic-ray spectrum can be properly included. The 
nuclei are injected from the maximum redshift z max = 6 
to the minimum redshift z m j n = 0.001 (corresponding to 
a distance D = 4 Mpc). The number density of hyper- 
novae at redshift z is assumed to follow the star forma- 
tion history (SFH), given by (Hopkins & Beacom 2006; 
Yiiksel et al. 2008) 

r (i + z) 3 -v< i 

n(z) cx I (1 + z)-°- 3 ,l < z < 4 (9) 
[ (l + z)- 3 - 5 ,z > 4. 

We collect all the nuclei and secondary products that 
arrive at the Earth (z = 0) and put them into the cor- 
responding energy bin and species bin according to their 
Lorentz factors j(z = 0) and mass numbers A(z = 0). 
Then we can get the energy spectrum of cosmic rays in 
each species bin, or get the all-particle energy spectrum 
by adding up all species particles at a certain energy. 

3.2. The final spectrum and composition 

The cosmic ray composition at the sources strongly 
affect the final composition and energy spectrum of cos- 
mic rays arriving at the Earth. First, we consider the 
chemical composition at the source equal to that of the 
WR stellar wind, where the mass ratio among dominant 
elements is M Ho : M c : M = 0.32 : 0.39 : 0.25 with neg- 
ligible H. The initial spectrum at the source is assumed 
to be a power law with p = 2. The result of the final spec- 
trum and composition are presented in Fig. 1. One can 

7 http://www.talys.edu/ 

8 Since the neutrons emitted by the photodisintcgration process 
will decay into protons very soon compared to the propagation 
time, we just treat the neutrons emitted as protons of the same 
energy in the calculation. 



see that with the typical WR stellar wind composition, 
the cosmic ray flux drops too fast at high energies to ac- 
count for the observations by PAO. The apparent reason 
for the fast drop is that the low E max for intermediate 
mass nuclei such as C and O results in an exponentially 
cutoff at energy < 10 20 eV. However, the situation will 
not be ameliorated even if we raise -E ma x,o to higher en- 
ergies, e.g. to 3 x 10 20 cV. This is because that lower en- 
ergy background photons, which are more abundant, are 
involved in the photodisintegration interaction for inter- 
mediate mass nuclei due to their higher Lorentz factors, 
compared with heavy nuclei of the same energy. There- 
fore the attenuation lengths for the intermediate nuclei 
with energies larger than 10 20 eV are so small (e.g. Puget 
et al. 1976; Allard et al. 2006) that even those who emit- 
ted by the nearest sources at z min = 0.001 (i.e. 4Mpc) 
will be effectively attenuated by the background pho- 
tons and hence the flux of cosmic rays at such energies is 
severely suppressed after propagation. Only if -B m a X ,o is 
extremely high (e.g. ~ 3 x 10 21 eV) so that secondary pro- 
tons of energies up to 2 — 3 x 10 20 eV can be produced in 
photodisintegration process, and then the spectrum be- 
comes flat at the highest energy end (Allard et al. 2007). 
However, such a high -E max ,o can hardly be achieved in 
the hypernova scenario. 

Presence of heavier nuclei in the source composition 
would be beneficial to fit the observation data, since 
heavy nuclei have higher steepening energies due to in- 
teractions with background photons. Hypernova ejecta 
produced after the explosion can provide heavy elements 
besides intermediate mass elements. Fig. 2 shows the 
propagated spectrum for cosmic ray composition at the 
sources equal to the hypernova ejecta composition in 
the "CO138E50" model for SN 1998bw (Nakamura et 
al. 2001). With the contribution by heavy nuclei at the 
sources, the propagated spectrum is in better agreement 
with the PAO data than the wind composition case. 

The fit can be further improved if some non-standard 
effects are taken into account. In panel a of Fig. 3, we 
increase the initial iron abundance by a factor of 3, i.e. 
iron nuclei constitute a fraction of ~ 25% of the total 
mass at the sources. This could be achieved since the en- 
richment of heavy elements is not uniformly distributed 
spatially, as already seen in some Galactic supernova 
remnants (e.g. Hwang et al. 2000). A larger explosion 
energy would also lead to more heavy elements synthe- 
sized during the explosion (e.g. Metzger et al. 2011). As 
iron nuclei have larger steepening energy in the spectrum 
due to interactions with background photons, a harder 
spectrum is expected at the highest energy if more iron 
nuclei are present. This is consistent with earlier results 
in Allard et al. (2008). Another possible effect is a hard 
injection CR spectrum. Although in conventional shock 
acceleration theory, the power-law index is p ~ 2 , some 
observations (e.g. Horns & Aharonian 2004; Green 2009; 
Reynolds et al. 2011; Abdo et al. 2011) and theoretical 
calculations (e.g. Ellison et al. 1996; Malkov 1997; Vainio 
& Schlickeiser 1999; Vainio et al. 2003; Tammi & Duffy 
2009) have suggested harder cosmic-ray spectra with in- 
dices p < 2. In panel b, we show the propagated spectrum 
for an injection spectrum with power-law index p = 1.6. 
One can sec that the fit to the data gets better than 
the case of p = 2. Besides, the source density distribu- 
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tion in the nearby universe could be nonuniform. Since 
our Galaxy locates inside the Local group and Local su- 
percluster, we examine the effect of local overdensity of 
sources on the propagated spectrum. Panel c shows the 
spectrum for a local overdensity in the source number 
by a factor of 2 relative to the average within the size 
of 30 Mpc (Blanton ct al. 2001). Due to the attenua- 
tion by background photons, higher energy cosmic rays 
are mainly contributed by closer sources. Thus a local 
excess of the source number density increases the flux of 
higher energy cosmic rays and thus hardens the propa- 
gated spectrum at the highest energy end. 

We should note that, although the deviation between 
the propagated spectrum and observational data is rela- 
tively large in the wind composition scenario, the situa- 
tion can also be improved if some of the above effects are 
taken into account. Panel d of Fig. 3 presents the result 
after considering a harder injection cosmic-ray spectrum 
of p = —1.6 as well as a local overdensity in the source 
number by a factor of 2. The theoretical energy spectrum 
agrees reasonably well with the observational data. 

In Fig. 4, we present how the average mass number 
of cosmic rays evolves with energy for different scenar- 
ios discussed above. Note that the starting energy of this 
figure is set to 3 x 10 18 eV since our model do not account 
for the CR spectrum data below this energy, for which 
a Galactic component contribution is needed. One can 
see that, for the hypernova ejecta composition scenarios, 
the average mass number increases gradually with en- 
ergy, which is consistent with the finding by PAO that 
the composition of UHECRs becomes increasingly heavy 
with energy. In the wind composition scenarios, the av- 
erage mass number increases more slowly with energy. 
However, since the hadronic interaction models at such 
high energies are not well understood, such a composi- 
tion is still consistent with the measurements within the 
uncertainties of theoretical expectations. 

4. DISCUSSIONS AND CONCLUSIONS 

We have shown that, with some fraction of heavy nu- 
clei such as iron in the source composition, the semi- 
relativistic hypernova model can explain the energy spec- 
trum and composition of ultra-high energy cosmic rays 
as measured by PAO. The heavy elements may origi- 
nate from the hypernova ejecta and are brought into the 
forward shock region via the Rayleigh- Taylor (R-T) in- 
stability. R-T fingers are seen to be able to reach the 
forward shock front or overtake the forward shocks (e.g. 
Jun et al. 1996; Blondin & Ellison 2001; Blondin et al. 
2001; Rakowski et al. 2011), so they should be able to 
bring heavy elements into the ambient medium. Obser- 
vations of some Galactic supernova remnants such as Cas 
A and SN 1006, as well as SN 1987A in Large Magel- 
lanic Cloud, have provided direct and indirect evidences 
for protrusion of heavy elements in the material beyond 
the blast wave (e.g. Erickson ct al. 1988; Fesen & Gun- 
derson 1996; Hwang et al. 1998; Hughes et al. 2000). 



When the ejecta is being decelerated by the swept- 
up ambient wind material, reverse shocks are also form- 
ing. The reverse shock may accelerate particles in the 
shocked ejecta, which provides an alternative site where 
UHE heavy nuclei are accelerated. The reverser shock in 
the hypernova case is typically mildly relativistic with a 
speed about /3 rs ~ 0.5, as shown in the appendix. Since 
the internal energy density in the shocked ejecta and 
shocked wind are equal, the magnetic field in the reverse 
shock is £> rs = (esr/es) 1,/2 -B, where es r is the magnetic 
field energy equipartition factor in reverse shock and €b 
is the same factor for the forward shock. Then we ob- 
tain the maximum energy of particles accelerated by the 
reverse shock 



= 1.8 x 10 



20 



1/2 



0.5 



(10) 



which can also reach ultra-high energies for heavy nuclei. 

In summary, hypernovae remnants are shown to be po- 
tential sources of UHECRs. They have been discovered 
to be present within the GZK horizon, ~ 100 Mpc. Due 
to the semi-relativistic ejecta produced by the explosion, 
they can accelerate particles to ultra-high energies. Hy- 
pernovae naturally provide intermediate mass or heavy 
elements either in the form of stellar wind or the hyper- 
nova ejecta itself. In this paper, we have calculated the 
propagated spectrum and composition of UHECRs arriv- 
ing at the Earth given an initial cosmic ray spectrum and 
composition at the sources produced by hyerpnovae. We 
find that with the source cosmic-ray composition equal 
to that of the hypernova ejecta, the propagated energy 
spectrum and composition are compatible with the mea- 
surements by PAO. 

There are suggestions that past Galactic hypernova 
could contribute to cosmic rays in the energy range above 
the knee (Budnik et al. 2008) or even UHECRs (Calvez 
et al. 2010). However, the event rate of such Galactic hy- 
pernovae is largely unknown. There has been evidence 
indicating that hypernovae associated with sub-energetic 
GRBs are preferentially found in low-metallicity galaxies 
(e.g. Woosley & Bloom 2006; Stanek et al. 2006), which 
would imply that hypernova rates in normal metallicity 
galaxies such as our Milky Way may be low. 
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APPENDIX 
THE SPEED OF REVERSE SHOCKS 

When the ejecta is being decelerated by the swept-up ambient wind material, reverse shocks are also forming. 
The shocked wind and shocked ejecta are in pressure balance and separated by the contact discontinuity. The semi- 
relativistic ejecta get decelerated at _Rhn — 5 x 10 16 £'fc,5i(ro/2) _2 Mr5 Wt«,3cm, where Ek and T are, respectively, the 
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total kinetic energy and the initial Lorentz factor of the semi-relativistic ejecta. Since the shock expansion is expected 
to be adiabatic, we have 

E s h,w + -E-sh,oj = Ek, (Al) 

where E s ^ w ~ 4:/3na(3 2 T 2 R 3 p w (R)c 2 (Blandford & McKee 1976) is the total energy transferred to the shocked wind 
from the blast wave, -E s h,ej — T 2 (e' + p'c 2 + /3 2 p')V^h,cj is the total energy in the shocked ejecta, T is the Lorentz 
factor of the shocked material and p w (R) is the wind density at radius R. Here a is a coefficient with a — > 0.35 when 
(3 — V 1 and a — \ 0.73 when /3 — > 0. We use <r = 0.73 — 0.38/3 to bridge the gap between the non-relativistic case and 
ultra-relativistic case, as used in Huang et al. (1998). e' = ^T"^ 1 (7 rs — l)p c jc 2 , p' = ^"^V ej an d p' = (7 — l) e ' are 
respectively, the internal energy density, rest mass density and pressure in the comoving frame of shocked material, 
and 7 is the adiabatic index. V s h, C j = Vo/( 7 ^1^" ) is the volume of shocked ejecta, where Vq is the volume of the ejecta 
before being shocked. Here 7 rs is the Lorentz factor of the reverse shock, which relates with the Lorentz factor L of 
the shocked material by 7 rs ~ r L(l — flofise), where /3q and /3 se are the velocities of the unshocked ejecta and shocked 
ejecta respectively. Solving Eq. (Al) with typical values for the parameters, we finally get the Lorentz factor and 
speed of the reverse shock, 

7«1.15,ft B = \/l - I/7 2 « 0.5. (A2) 
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Fig. 1. — The propagated spectrum of cosmic rays assuming a source composition equal to that of the stellar wind of Wolf- 
Rayet stars, where the mass ratios of dominant elements are He: C: 0=0.32: 0.39: 0.25, and a flat source spectrum with 
p = 2 (see text for more details). The black solid line represents the all- particle flux, while other lines represent the contribu- 
tions by different species as shown in the legend. The black squares are the most recent data measured by PAO, taken from 
http: / /www. auger.org/technical_info/ICRC2011/CombincdSpectrum_Auger2011.txt. 
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Fig. 3. — Panels a, b, c show the propagated cosmic ray spectrum for the source composition equal to that of the hypernova ejecta, similar 
to Fig. 2, but with some modifications taken into account. In panel a, we increase the abundance of iron nuclei in the ejecta by a factor of 
3. In panel b, we set a hard initial spectrum with p = 1.6. In panel c, we assume that the local source number density is two times higher 
than the average within a distance of 30 Mpc. The dashed line in panel c represents the contribution by sources beyond 30 Mpc while the 
dash-dotted line represents the contribution by local sources within 30 Mpc. Panel d shows the propagated cosmic ray spectrum for the 
source composition equal to that of stellar wind, but with the effects of a hard spectrum and a local overdensity being taken into account. 
See text for more details. 
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Fig. 4. — The average mass number of propagated cosmic rays arriving at the Earth as a function of the cosmic-ray energy for different 
cases shown in the above figures. . 



